Bounded Global Optimization for Polynomial Programming 
using Binary Reformulation and Linearization 



Joseph W. Norman 
University of Michigan 
jwnormanOumich. edu 

May, 2012 
ABSTRACT 

This paper describes an approximate method for global optimization of polynomial programming 
problems with bounded variables. The method uses a reformulation and linearization technique to 
transform the original polynomial optimization problem into a pair of mixed binary-linear programs. The 
solutions to these two integer-linear reformulations provide upper and lower bounds on the global solution 
to the original polynomial program. The tightness of these bounds, the error in approximating each 
polynomial expression, and the number of constraints that must be added in the process of reformulation 
all depend on the error tolerance specified by the user for each variable in the original polynomial 
program. As these error tolerances approach zero the size of the reformulated programs increases and 
the calculated interval bounds converge to the true global solution. 

1 INTRODUCTION 

This paper describes a method for approximating with interval bounds the global optimum of a polynomial 
program with bounded real variables and no special convexity properties. Such problems arise in many 
contexts including the analysis of parametric probability models [8] H] . For that application and many others 
it is sufficient to generate interval bounds on the global minimum and maximum solutions of the polynomial 
optimization problems involved; it is not required to compute those solutions with infinite precision. The 
successive approximation method presented here uses a reformulation and linearization technique to create 
a hierarchy of pairs of mixed binary linear programs whose solutions bound the true global optimum of the 
original polynomial problem. 

The method allows the user to adjust directly and incrementally how many variables and constraints 
are added in the reformulation stage, and to compute hard bounds on the error in approximating each 
polynomial before solving the reformulated problem; thus the user can control computation to trade the 
tightness of the generated interval bounds against the time required to compute them. The reformulated 
problems can be solved by well-developed branch-and-bound or branch-and-cut methods for mixed integer 
linear programming for which codes are available in commercial and free optimization software. The general 
global polynomial optimization problem is NP-hard and it remains an open question how well the bounds 
provided by the proposed method converge in practice and how best to control computation to produce the 
most useful answers with the smallest amount of computation. 

The main contributions of this method are control of computation and definitive interval-based answers 
with predictable and explicit bounds on potential errors. Based on the parameters provided by the user, 
the method computes bounds on the possible error in the linear approximation of the objective and each 
constraint, before the reformulated optimization problem is solved. This predictability of error bounds means 
that the method can generate hard interval bounds on the global solution to each polynomial optimization 
problem, conditioned on the feasibility of that problem. Feasibility can sometimes be confirmed or refuted by 
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solving the reformulated problems; and even when feasibility cannot be established definitively, the method 
offers useful linear bounds on the set of points that are potentially feasible. 

1.1 Preliminaries 

Let us consider a polynomial program PP: 
minimize f (xi , . . . , x n ) 

subject to gi (xi , . . . , x n ) 0, g2(xi, • • • ,x n ) < 0, . . . , g q (xi, . . . ,x n ) < (1) 
and oci ^ X] ^ |3i , a.2 ^5 y-i ^ p2, • • • , <x n ^ x n ^ |3 n 

in which the objective f : R n — > R is a polynomial function of the variables Xi , . . . ,x n as is each constraint 
g,. The lower bound on each real variable xt is at and the upper bound is |3i. Each variable Xi may take 
positive, negative, or zero values. Fractional objectives can be accommodated using the Charnes-Cooper 
transformation [T] . It will be demonstrated how to generate from PP a pair of mixed integer- linear programs 
such that the global minimum of PP is bounded by the solution to each linear approximation. 

Let us use x = (xi , . . . ,x n ) to represent the tuple of variables used in the polynomial program PP. Let 
us use O x to represent the set of values of x that satisfy the bound constraints included in the polynomial 
program PP: 

fx = {(Xl, . . .,X n ) : OCI < Xl < Pl,...,On < Xn < PtJ (2) 

The bound constraints O and the polynomial constraints gj ^ define a semialgebraic set of feasible values 
for x. Thus we can rewrite the polynomial optimization problem PP as follows: 

minimize f(x) 

subject to gi (x) < 0, g2(x) ^ 0, . . . , g q (x) ^ (3) 
and x e <D X 

Let us use fpp(x*) to denote the true global minimum solution to the problem PP, where x* identifies a 
feasible (though not necessarily unique) point at which that minimum occurs. The goal is to develop a pair 
of mixed integer-linear programs LP~ and LP + using new variables w with the property that their solutions 
bound the true global minimum: 

f LP -(w*) ^ f PP (x*) s; f LP+ (w*) (4) 

It is also desirable for the linear approximations to provide trustworthy information about the feasibility 
of the original nonlinear program, reporting whether feasibility has been confirmed, refuted, or not yet 
determined. The reformulation and linearization method presented below accomplishes these goals. 

The tuple w of variables used in the reformulated programs contains three classes of variables: 'unit 
variables' denoted ui,U2, . . . ,14, of which each binary variable 6 {0,1} is allowed to take the integer 
value zero or one; 'remainder variables' denoted ri,T2, ... r n of which each real variable ^ T\. ^ 1; and 
'unit-product variables' denoted yi ,ij2> • • • i^\\> of which each real variable < yi ^ 1 . The number of unit 
variables is <t> and the number of unit-product variables il>. Therefore the total number of variables in w 
used in reformulated programs is n + 4> +il>, where n is the number of variables in the original program PP. 
Thus the reformulated variables w are given by: 

w = (u 1 ,u 2 ,...,u (t) ;r 1 ,r2,...,r n ;u 1 ,'y2,.--,'y^) (5) 

As part of the reformulation process some number p of linear constraints are constructed and added to the 
problem; these constraints define a feasible set <P W of values for the reformulated variables in w. Thus the 
problem is lifted from the space M n to the space M n + t t>+ 1 ^ and linearized in that space. Analogous to the 
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bounding box <D X that limits the feasible values of x € R n a new bounding box O w limits the feasible values 

of w e R n+ * + ^. 

In order to control the linearization and reformulation process the user supplies parameters <Ti , csz, ■ ■ ■ , d n 
where each o~i determines how the corresponding original variable Xi is to be approximated. The tightness 
of the interval bounds computed on the global solution to PP, as well as the number ct> of unit variables, the 
number i|) of unit-product variables, and the number p of constraints that must be used in the reformulation 
process, all depend on these user-supplied parameters. 

2 PRIOR WORK 

2.1 Reformulation and Linearization 

The method presented here is based on the principles set out by Li and Chang [7 , which in turn reference the 
reformulation linearization technique of Sherali and Tuncbilek [11 ; it is also a progression of the linearization 
and bounding techniques described in the author's doctoral dissertation [8j. The main extension relative to 
the work of Li and Chang is that here there is a defined relationship between the solution to each linear 
reformulation and the true solution to the original polynomial program: the minimum solution of one linear 
program is less than or equal to the true global minimum, which is in turn less than or equal to the solution 
of the other linear program. The original work presented in |7J offers the exact solution to a problem related 
to the original polynomial problem, but with the exact qualifications on that relationship left unspecified. 
These desired solution properties are produced by explicit calculation of the possible error in approximating 
each polynomial expression through linearization combined with reasoning about the inequalities involved. 
Additionally, in this work a different method is used to linearize products of two or more variables. Finally, 
the number of variables and constraints necessary to reformulate an optimization problem by this method 
is calculated. As in the original work no assumptions are made about the convexity of the feasible set or 
objective function in the polynomial problem; hence the technique presented here works for maximization as 
well as minimization of objectives. For simplicity of presentation only minimization problems are described 
below. 

The reformulation-linearization technique (RLT) described by Sherali and Tuncbilek [TT] presents an 
approach to solving general polynomial programming problems which is quite similar. In each case the 
problem is lifted from its original n-dimensional space R n to a higher-dimensional space via the addition of 
new variables, and linear constraints are imposed in that higher-dimensional space to provide approximations; 
the original RLT provides only outer bounds (e.g. upper bounds on the global maximum), whereas the 
present work provides inner bounds as well. Mechanically, the present work uses two kinds of linearization 
constraints (let us call them real-product-mean and binary-product-sum) which are different from the 
bound-factor, constraint- factor, and convex-variable-bounding constraints used in the traditional RLT. 
In Sherali and Tuncbilek's work the new variables in the lifted space map one-to-one to the monomial 
terms that comprise the basis of the original polynomial space; thus the number of new variables is fixed. 
Depending on which classes of linearization constraints are added, the number of additional constraints 
varies; the user may choose to add whole classes of constraints or select one at a time in an ad-hoc way. 
However in the present work the number of additional variables can vary as well as the number of additional 
constraints. Furthermore, these are both determined in a strict algorithmic way according to a small set of 
parameters provided by the user. 

One benefit of the approach taken here is that the user can make a sequence of reformulations at different 
granularities or resolutions. Moreover, the reformulation and subsequent linearization are separated from 
the branch-and-bound search used to solve the resulting linear programs (which are mixed integer linear 
programs in my case); that is, the branch/bound requirements are encoded in integer variables. This allows 
the current method to take advantage of well-developed codes for solving mixed integer linear programs. 
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This also allows the user to control the size of each reformulated program and thus the time required to 
solve it. We shall present in more detail the relation to the traditional reformulation linearization technique. 

2.2 Semidefinite Programming and Real Algebraic Geometry 

Besides the linear relaxations of the RLT, another option is to use semidefinite program relaxations of the 
original polynomial program, to compute a lower bound on the true minimum. As described by Lasserre 
[5] [6] this approach is in fact a generalization of the linear relaxation approach and is related to it through 
the theory of moments and its dual theory of representation of polynomials which have positive values over 
a semialgebraic set. Additionally, Floudas and Visweswaran [2] [3] present a technique to transform bilinear, 
quadratic, or polynomial programs to a new problem with partitioned variables and a certain convexity 
property; the transformed problem can be solved by a primal and relaxed dual approach. 

There are many methods to solve special cases of the polynomial programming problem, e.g. those in 
which all polynomial functions have positive values or other instances in which favorable convexity conditions 
are met. For the broader class of constrained global optimization problems there are also a variety of methods, 
including outer approximation and branch-and-bound techniques, as reviewed by Horst and Tuy [4 . In the 
classification of Horst and Tuy, the present method is a successive approximation method. It is deterministic 
rather than stochastic. The inner or outer bounds at each stage are computed by relaxation of the original 
polynomial program to a mixed binary linear program. Branching and bounding can be carried out to solve 
each of these linear programs; however the branching and bounding occurs in the (binary part of) lifted 
space not the original space of the polynomial problem. The current method uses only the primal form of 
the polynomial program and does not incorporate the dual problem in its analysis. 

Most successive approximation methods reformulate the problem in one step into a convex (linear or 
semidefinite) form; here we use a nonconvex intermediate (a mixed 0-1 integer linear program) which is 
NP-hard to solve but for which a great deal of work has been done to make algorithms that perform well on 
average. Note that successive approximation with semidefinite relaxations is one of the possible approaches 
to solving mixed binary linear programs which is distinct from the more common branch-and-bound or 
branch-and-cut methods. 

It is not clear to me if one can recover from each SDP relaxation any bound on the error in the approxi- 
mation (perhaps using the primal and dual solutions provides appropriate bounds on the global minimum), 
or if it is possible to determine error bounds on the approximation before solving it; it seems instead that 
after solving a particular relaxation it may be verified that it is indeed a global solution. 

3 BINARY REFORMULATION AND LINEARIZATION 

Here we develop the method to create a pair of mixed binary linear programs to bound the global solution 
to a polynomial program. 

3.1 Binary Expansion and Basic Inequalities 

Proposition 1 (Reformulation with unit variables) Consider a real-valued variable x t bounded 
by the constants oci ^ Xi ^ (3i. Given a positive constant Ki satisfying < $C (3t — oci it is possible 
to represent X| as a sum that involves some number <Ji of binary variables Uij € {0,1} and a 
nonnegative real variable rt € [0, 1]; 




(6) 
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Let us call each u^j a unit variable and v\ the remainder variable; each Ki is the error limit on the 
corresponding variable xi. In general it is necessary to add the constraint that the sum satisfies the 
original upper bound pi.: 

0d + iCi222l- 1 uy + iCiT i < Pi (?) 

i=i 

Otherwise it may be possible to choose values of Ki, uij and ri that would produce a value greater 
than (3. 

This first proposition is a restatement of Equation 2.1 in [7J, with additional detail provided here. The 
number o"i of unit variables required to reformulate a variable xi in this way is related to the error limit Ki.: 



i i Pi Wi c 

log 2 h 6i 

Ki 



(8) 



where the brackets indicate rounding up and Si = 1 if Xi is discrete or if Xi is continuous. The lowest 
possible error limit k* for a given number o"i of unit variables is given by: 

< - IBt < 9 > 

where again 8i = 1 if Xi is discrete. For each variable Xi the user may specify the number o"i of unit variables 
and have the system compute the corresponding error limit k* using Equation[9l or the user may specify the 
desired error limit Ki. and have the system compute the number of unit variables required by Equation [51 

Note that in the case that o"i = the error limit Ki = pi — <Xi and no unit variables will be added; 
note also that a discrete variable with cri = must have 0L\ — Pi in which case Ki is not needed for 
reformulation. Furthermore the error limit Ki can be zero only for a fixed variable x^ with bounds oci = |3i. 
Also, it turns out that it is not necessary to add any unit variables for a variable Xi which appears only 
linearly in PP; in this case we will take o"i = but Ki = pi — oci (thus the reformulation of such an Xi is 
KiTi which is equal to (Pi — <X\]Ti with ^ T\ ^ 1). Finally, we can model a discrete variable with steps 
Ki using the substitution above but omitting the remainder variable Tij such a variable would take values 

{(Xi, OCi + Ki, Oti + 2k u . . . , OCi + CTiKi} with OCi + OiKi < Pi. 

Example 1 Consider the variables x = (xi,X2,X3) bounded by 2 ^ xi ^ 5, ^ X2 ^ 10, and 4 ^ X3 ^ 8. 
Let us use three unit variables to reformulate xi and two unit variables to reformulate each of X2 and X3. 
In other words o"i = 3, 02 = 2, and ff3 = 2. According to Equation [9] the smallest possible error limits are 
then Kf = 0.375, k 2 = 2.5, and K3 = 1. The reformulated variables according to Proposition [l] are: 

X! = 2 + 0.375ui ,i + 0.75u! >2 + 1 .5ui ,3 + 0.375r! 

x 2 = 2.5u 2 ,i +5u 2 ,2 +2.5r 2 (10) 
x 3 = 4 + u 3>1 +2u 3i2 +r 3 

With the following constraints added by Equation [7J 

2 + 0.375x11,1 +0.75x11,2 + 1.5ui >3 + 0.375T! < 5 

2.5u 2) i +5u2,2+2.5r 2 < 10 (11) 
4 + u 3 ,i +2u 3 , 2 +r 3 sC 8 

Proposition 2 (Linear bounds on real products) Consider nonnegative real-valued variables ri,r '2, . 
each bounded by ^ T-t ^ 1/ assume n > 0. The difference between the mean and the product of these 
bounded variables is limited by the number of variables: 

< ( ri+r2+ n - +T j -(i,r 2 .., n ) < (^1) (12) 



■ • > ' n 
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which can be rewritten as: 



< fa +T 2 + --- + T n )-n(T 1 r 2 »-T n ) < n-1 (13) 

This proposition is a generalization of Equation 2.5 in with the following proof added here. It is clear 
that the relationships are true when every rt is zero, in which case both the mean and product are zero; and 
when every r\ is one, in which case both the mean and product are one. For values other than these the 
greatest difference between the mean and the product will occur when one of the variables is zero and the 
others are one; in that case the mean will be (n— 1 )/n and the product will be zero, satisfying the equations 
above. 

Proposition 3 (Linear bounds on binary products) Consider binary variables ui ,U2, . . . ,u n with 
each ui € {0, 1} and also one real-valued 'remainder' variable r that satisfies ^ r ^ 1. The follow- 
ing relationships hold between the product of the binary variables, the remainder variable, and the 
individual binary variables: 

. , u,u 2 ---u n r ^ r + U!+u 2 H +u n -n 

uiu 2 ---u n r ^ m, x = 1,...,n (14) 

ui u 2 • • • u n r ^ r 

If the real variable r is omitted the appropriate relationships between the sum and product of the unit 
variables are: 

. . -i U!U 2 • • -Un > 1 +ui +u 2 H hu„-n 

UiU 2 ---Un s$ i=1,...,u . , (15) 

Ul U 2 • • • Un < 1 

This third proposition is a restatement of Proposition 1 from [7J and the proof that appears there. It is clear 
that if any unit variable = then the product ui u 2 • • • u n = 0. In that case the sum ui +u 2 + • • • +u n ^ 
n — 1 as one of the terms is zero and each of the others is not greater than one. The constraints above are 
satisfied in this case. If every unit variable uj = 1 then the sum ui + u 2 + • • • + u n = n and the product 
uiu 2 • • -u n = 1, in which case the constraints are also satisfied. 

Let us introduce for this discussion a unit-product variable y — ui u 2 • • • u n r to represent the product of 
the binary variables ut through u n and the real- valued variable r described above; we can use y = ui u 2 • • • u n 
for the product of the binary variables alone when no real variable r is included. Each unit-product variable 
y is continuous and bounded by ^ y ^ 1 . 



3.2 Reformulation and Linearization of Polynomials 

Now we have the tools to reformulate any polynomial expression from its native x-variables into a linear 
function of the unit variables u, the unit-product variables y, and the remainder variables r. This procedure 
is a generalization of Propositions 2 and 3 from |7J which describe only the linearization of products of two 
or three variables. 

3.2.1 General Polynomial Form 

Consider that any polynomial function g(xi , x 2 , . . . , x n ) can be represented as the sum of several terms 
Ckrak where each term is the product of a real coefficient Ck and a monomial m.k, the latter of which is the 
product of several x-variables: 

g(xi,x 2 , .. .,x n ) = Cimi+c 2 m 2 H hc t m t , c k 6l, ra k = Y\ Xi ( 16 ) 

iei k 

*There is a typographical error in Equation 2.5 in [7|: the final term should be jU)]U)2 instead of lu>ju>2. 
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Let us adopt the convention that the first monomial mi = 1 so that the coefficient Cj is the constant term in 
the polynomial g. Each tuple Ik of indices identifies the original x-variables included in the product for the 
monomial m.k- An index i may occur in Ik more than once; if a tuple Ik is empty then the corresponding 
product of zero variables is taken to be unity. Note that the size |Ik| of the tuple of indices for a monomial 
mic is the degree of that monomial. Let us use m to represent the list of all monomials which occur in a 
particular optimization problem PP, and t for the number of such monomials. 

Example 2 Consider the following polynomial program PP1 which is the example problem PP(O) presented 
in Sherali and Tuncbilek [llj and reproduced as Example 1 in Li and Chang [7|: 

minimize : 5x2 + *3 + — 2xi xz — 3xi X3 + 5x2X3 

-X3 +X!X 2 X 3 

subject to : 4xi + 3x2 + X3 ^ 20 
X] + 2x 2 + x 3 > 1 
and : 2 ^ xi ^ 5 
< x 2 10 

4 ^ x 3 s$ 8 

The list of monomial terms (basis) used in PP1 is: 

m = (l,X 1 ,X2,XiX2,X3,XiX2X3,xf,X 1 X3,X2X3,x|) (17) 

The first term mi = 1 has the empty tuple I] = () of indices and degree 0; the last term m.10 = x\ has the 
tuple 1 1 o = (3,3) of indices and degree 2. The number of terms in m is t = 10 



3.2.2 Reformulation of Monomials with Sums of Unit Variables 

Now consider just one monomial rat = Oiei k x i °f a polynomial g as given in Equation fl6l Using Propo- 
sition [l] we can substitute a sum involving several unit variables and a remainder variable for each original 
variable Xi included in the product that defines this monomial m.k where each at is the lower bound of the 
corresponding variable Xi and each o~i is the number of unit variables necessary to represent xt within a 
tolerance of iq,. Carrying out the multiplication to distribute the product over the sum yields an expression 
for the monomial m.k of the form: 

mk = ]^[ ai+K i ^2 H1 u i) j + K i r i = z k) i+Zk,2H r-Zk,s k (18) 

ieik \ j=i / 

where each Zk,t designates the i-th element in the sum that constitutes the distributed product. With 
reference to Equation [18] the number Sk of such elements is limited by: 

as the sum used to replace each variable Xi according to Proposition [l] contains o"i unit variables, one 
remainder variable, and one constant lower bound <Xi if that bound is nonzero. 

Now, it is clear from the structure of the inner sum and outer product shown in Equation IT8l that every 
element Zk,£ in the sum must be the product of some numbers of constant lower bounds a, unit variables 
u, and remainder variables r: 

Zk,« = (aijaq • ■•a^J (Ki^'-V,,}, Ki 2 2 ]2 -V i2)j2 • • • K inu 2'"-- 1 u inu)jnu ) (K h ,r h , K h2 r h2 • • • k^t^O) 
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where n a is the number of constants, n u the number of unit variables, and n r the number of remainder 
variables. Note with reference to Equation [18] that the number of each kind of element is limited by the 
degree die = |Ikl of the monomial mt: 

n„ ^ d k , n u ^ d k , n r ^ d k (21) 
Let us introduce the constant a.k,e to simplify notation: 

= ^(ai;^---ai;J(K il K i2 ...K inu 2^- 1 2^- 1 ..-2i^- 1 )(K hl K h2 ...K hnT ) (22) 

where we take the product of an empty set of items to be unity and where we omit division by n r if n r = 0. 
With this constant a^e we can rewrite the element Zk,£ shown in Equation [20l as: 

= a^Tir (u iuh Ui 2th •••u illu)jnu ) (r hl TH 2 ■■■r KnT ) (23) 

where n r is omitted if it is zero. Note that cik,£ may be positive or negative. 

Example 3 The ninth term in the list m of monomial terms used in the problem PP1 shown in EquationfTTl 
is vn<? — x 2 x 3 . Substituting the reformulation for each variable given in Equation [TO] this term becomes: 

m.9 = (2.5u 2 ,i + 5u 2 ,2 + 2.5r 2 ) (4 + u 3 j + 2u 3)2 + r 3 ) (24) 

Carrying out the multiplication gives the following sum of 12 terms: 

m.9 = 10U2,1 +2.5U2,1U3 ) 1 +5U2,lU3 i 2 +2.5U2,lT3 +20U2,2 +5U2,2U3 ) i (25) 
+ 10U2,2U 3 ,2 +5U2,2l"3 + 1 0X2 +2.5u 3 ,lT2 + 5u 3 ,2T2 + 2.5r2r 3 

These terms are the elements Z9,i through Z9,i 2 as described in Equation [TU Let us examine a few of these 
elements. For the fourth element z?^ — 2.5u2ju 3 j we have: the number of unit variables n u = 2, the 
number of remainder variables n r = 0, and the constant 0.9,2 = 2.5. For the sixth element Z9,6 = 2.5u2,ir 3 
we have: n u = 1, n r = 1, and a<?^ — 2.5. For the last element Z9 12 = 2.5r 2 r 3 we have: n u = 0, n r = 2, 
and a 9 ,i2 = j{2.5) = 1.25. 

3.2.3 Linearization of Products of Remainder and Unit Variables 

We can use the propositions above to linearize each element Zk,£ of the sum that constitutes the reformulated 
monomial term mk in the following way. First let us introduce the notation [zk,« J with double square brackets 
to denote the linearized version of Zk,£ • Let us begin to compute this linearized form by replacing the product 
of remainder variables in Zk,f with their mean: 

t t 1 ^ /r Kl +r h2 + --- + r hn \ 

[zic.ij = aic )£ n T (utuhUi^ •••UWu.jnJ I J ( 26 ) 

Simplifying this expression yields: 

[z k ,i] = a k){ (ui,,j,Ui 2 ,j 2 ■■■u iriutUu ) (r Hl +r h2 +---+Th Tlr ) (27) 

Note that if the element Zk,{ contains no remainder variables (n r = 0) the linearized form is unchanged 
from the original. Let us now replace the product of the unit variables (u-l, ,j , u-i 2i j 2 • • , 'Ud T1 . u ,j TVu .) an d each 
remainder variable in Equation [27] with a new unit-product variable y k { h such that: 

yk,£,h = (ui^^Ui,,), •••u inuiKu )r h , h= 1,2,...,n T (28) 
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The linearized element Jzk,«] can now be written as weighted sum of several unit-product variables: 

[zk.t] = dk,« ("yk,«,i + yk,«,2 H (-y k>£jTtr ), ifn r >0 (29) 

In the case that there are no remainder variables in the element z k t an d thus n r = let us introduce a 
single unit-product variable: 

yk,<!,i = (•Ud,,j,ui 2 ,j 2 ---Ud nu) j nu ) ( 30 ) 

and express the linearized element [Zk,«] appropriately: 

[z k ,{] = a M y k , £ ,i, ifn T =0 (31) 

To simplify notation let us introduce the variable rik i which equals n r if n r > 0, or one otherwise. Thus 
each linearized element [z k ,{] can be expressed: 

it ti f 1 if Tu = 

[z k ,«l = a k)£ 2_ yic,«,H, n k)£ = . f > Q (32) 

h=1 ^ 

Example 4 Let us return to the elements that comprise the expanded form of the monomial term m.9 = X2X3 
described in Example [3] above. For the element z?,2 = 2.5u.2,iU3,i , the linearized form [29,2] = 2.5y9,2,i 
using the solitary unit-product variable y9,2,i = 112 ) 1 U3 1 1 introduced according to Equation [30l For the 
element Z9,4 = 2.5u2,iT3 the linearized form [2:9,4]] = 2.5y9,4,i using the unit-product variable 1)9,4,1 = 
u 2,iT"3- For the element Z9,i2 = 2.5t2T3 the constant (19,12 = 1-25 according to Equation [22] and the 
linearized form [Z9J2I = 1.25(r2 + Ts) according to Equation [27] According to Equation [28] we could 
introduce two trivial unit-product variables 1)9,12,1 = i"2 an d 1(9,12,2 = T3 and express the linearized term 
[^9,12] = 1-25 (1)9,12,1 +1)9,12,2) using the standard form given in Equation [32] 



3.2.4 Linear Constraints on Unit-Product Variables 

For each element z k ,{ of the sum shown in Equation [T8l one unit-product variable y k ,£,h. has been introduced 
for each of the n r remainder variables in Equation [27] (or a single product variable y k ,{,i if n r = 0). In order 
to satisfy Proposition [3] it is necessary to add the following constraints for each unit-product variable i) k ,£,H 
thus introduced: 

l)k,f,h ^ u ii,ji 
!)k,£,h ^ u i2,j2 

where Th refers to the remainder variable associated with the unit-product variable y k) { h and with the 
substitution of unity in the place of Th if n T = 0. The number of added constraints for each element z k ,e is 
tLt(t!u + !)• For a unit-product variable yk,t,H = Ti which equals some remainder variable (i.e. n r = 1 and 
n + u = 0) or a unit-product variable y k ,£,h = Ud,j which equals some single unit variable (i.e. n r = and 
n u = 1)i the constraints implied by Equation [33] are trivial and need not be added; it is sufficient to note 
the identities. In fact, in such cases the original variables can be used and new unit-product variables need 
not be introduced. 

Let us use O w to designate the set of values of w that satisfy the linear constraints shown in Equation [33l 
for all unit-product variables introduced, as well as the upper bound constraints as shown in Equation [7] for 
all reformulated original variables. Let us use p to designate the number constraints used to specify the set 



yk,£,h > Th+Ui,,], +U i2)j2 +"-+U init ,j nu 



(33) 
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Example 5 Let us use the elements discussed in Examples [3] and [4] to illustrate the constraints on unit- 
product variables. The complete linearization of the reformulated monomial term m.9 is given by: 

[m 9 ] = 10u 2 ,i +20u 2 , 2 + H.25r 2 + 1.25r 3 +2.5y 9 ,2,i +5y 9 ,3,i (34) 

+ 2.51)9,4,1 +51)9,6,! +10lj9,7,l + 5y 9)8 ,1 +2.5y 9 , 10 ,1 +51/9,11,1 

Compare this with the reformulated but not yet linearized term shown in Equation [25] note that the 
corresponding terms appear in different orders in the two equations. The linearized form above uses the 
following non-trivial unit-product variables: 



^9,2,1 


= U2,lU 3) l 


¥9,3,1 


= U 2 ,lU 3) 2 


y9,4,i 


= u 2 ,ir 3 


¥9,6,1 


= U2,2U 3 ,! 


¥9,7,1 


= U 2 ,2U 3 , 2 


¥9,8,1 


= U-2,2r 3 


¥9,10,1 


= u 3) ir 2 


¥9,11,1 


= u 3 , 2 r 2 



(35) 



These unit-product variables are accompanied by the constraints described in Equation [33] For example, 
for the unit-product variable 1)9,2,1 = 1^2,1^3,1 , Equation l33l implies the following constraints: 

u 2 ,i > 1)9,2,1 

u 3 ,i > 1)9,2,1 (36) 

U-2,1 +U 3 ,1 < 1 +1)9,2,1 

Similarly for the unit-product variable 1)9,4,1 = u 2 ,ir 3 Equation [331 implies: 

U2,1 1)9,4,1 

u 2 ,i +r 3 < 1 +1)9,4,1 (37) 

T"3 > ^9,4,1 

For 29,12 the constraints from Equation [33] on the unit-product variables 1^9,1 2,1 = 1*2 and 1)9,12,2 = t 3 are 
trivial and need not be added. 

3.2.5 Reformulation and Linearization of Entire Polynomials 

Recall from Equation [181 that a polynomial term m k representing the product of several variables xi can be 
represented as the sum of several elements z k ,{ 

Ta k = z k ,i+z k , 2 H hzic iSt (38) 

where each element z k ,£ is the product of some number of constants, unit variables, and remainder variables 
as shown in Equation [20] Using the reformulation technique above a linear approximation [z k ,{] can be 
generated for each element z k ,{ in this sum; adding these approximations yields a linear approximation 
[mic] for the polynomial term m k : 

[trtJ = [z k ,i] + [z k , 2 ] + --- + [z k , Sk ] (39) 

From the above and Equation [32] it is clear that the linear approximation for each monomial term m k is a 
weighted sum of unit-product variables ij k ,£,h.: 

Sk Tt k ,{ 

a k ,fu k ,«,h (40) 

£ = 1 H=1 
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where Sk is the number of elements in the sum and each rik^ is the number of remainder variables included 
in the product that defines each element, or 1 if there are none. The reformulation also requires constraints 
involving the unit-product variables yk,«,Ki the unit variables utj, and the remainder variables as shown 
in Equation [33] It is clear that this reformulation technique could be applied successively to each monomial 
CkTrik in a polynomial g as given in Equation [T6] to yield a linear approximation [g] of that polynomial: 

t s k n kjt 

[g] = ci ■ [mil + c 2 • I m 2l H h c t ■ [m t ] = XXX c * a *>*yM,H ( 41 ) 

k=1 «=1 h=l 

3.2.6 Bounds on Error from Linearization 

In this section we shall consider the error introduced by the linear reformulation process described above. 
The error is a function of the error limits Ki , k 2 , • • • > K n set by the user for the variables X\ used in the 
polynomial expression g. It follows from Equations [23] and [27] that for n r > 1 the difference between a 
linearized element [zk,e] and its true value Zk,e is exactly: 

[z M l]-z M = a k , { (u il0 -,u i2)j2 •••u illu)jn J ((r h , +r h , + ••• + r Kir ) -n r (r h ,r h2 •••r hn J) (42) 

Note that the constant dk,£ contains the product of several error limits K| as shown in Equation [22] In the 
case that n T = or n r = 1 no error is introduced; error is introduced only by linearizing elements with 
n r 2 that contain the products of two or more continuous remainder variables. 

According to Proposition [2] the difference [zk,{]] — Zk,£ shown in Equation [42] is bounded by: 

< [z k ,e] -z k ,f < a M( u i-i>ji u i-2,j2'" u iTv„,jTv u ) (^r — if a k ,{ > (43) 

As the product ("U-i, j, Ui 2) j 2 • • ■ u i rilt ,jtiu) °f the un ^ variables must be zero or one these error bounds are 
simplified to: 

s? [zk, £ ]]-Zk,e < a k ,f(n T -l), if a k ,f > 0, n r > (44) 
with the corresponding bounds for a negative coefficient Qk,£: 

aic,{(n r -l) ^ [zm! -z M ^ °> if ik^e < 0, n r > (45) 

Using these relationships it is possible to create a pair of functions to bound the error in the linear approx- 
imation of a polynomial. 

Proposition 4 (Scalar bounds on error in linear approximation) For a polynomial g as given 
in Equation ] lb] let us define the lower error bound E~(g) to be the amount by which the linearization 
[gj might underestimate the true value of g and the upper error bound E + (g) to be the amount by 
which the linearization [g] might overestimate the true value of g : 

E"(g] < [gj-g < E + (gl (46) 

With reference to the error bounds on elements Zk,£ given in Equations \44\ and \45\ the sum of unit- 
product variables yic,t,h that defines each linearized item [Zk«] shown in Equation\32\ and the sum of 
unit-product variables that defines each linearized polynomial [g] shown in Equation \41[ these error 
bounds E~(g) and E + (g) can be computed as follows. The lower error bound E~(g) is given by: 

E~(g) = I Y_ Xcka k ,£ (n M - 1) : c k a k ,f < i (47) 

lk=1 «=1 J 
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where the sum includes an element only when the product of coefficients Ck<ik,£ is negative. The 
corresponding expression for the upper error bound E + (g) is: 

H + (g) = |^^c k a M (n M -1) : c k a M >oj (48) 

lk=1 «=1 J 

including elements only for positive products of coefficients Ck<ik,£ . Note that the product Ck<ik,£ (n-k,£ — 1 ) 
is nonzero only when n^t > 1 (equivalently when n r > 1 ). 

Example 6 Returning to the element Z9J2 given in Example^ you can see from the reformulations for the 
individual variables X2 and X3 given in Example [T] that the element Z9J2 = 2.5r2i"3 is in fact the product 
( K 2i"2) (K3T3). Thus according to Equation [22] the constant — 2~ K 2 K 3- Equation [44] says that the 

difference between the linearized element and its original reformulation must be bounded by: 

< [Z9,12] -Z 9 ,12 < -jK2K 3 (2-l) 

Substituting K2 = 2.5 and K3 = 1 gives ^ [29,12] — 29,12 ^ 1.25. Error bounds for complete polynomials 
can be computed using Equations l47l and [481 For example the bounds on the objective function: 

f(xi,X2,x 3 ) = 5x 2 + x 3 + x] - 2xi x 2 - 3xi x 3 + 5x 2 x 3 

-%l +X!X 2 X 3 

of the problem PP1 in Example [2] turn out to be E"(f) = and E+(f) = 1.25. 



3.2.7 Linear Bounds on Polynomials 

The reformulation and linearization procedure above, along with the computed error bounds, allow us to 
compute linear expressions that provide upper and lower bounds on any polynomial g(x). These bounds 
are valid for all values of x within the feasible set O x . 

Proposition 5 (Linear bounds on polynomials) Consider a polynomial g which is a function of the 
variables x = (xi ,X2, . . . , x n ). Let O x denote the set of values of x that satisfy the bound constraints 
<X-i ^ xi ^ Pi as shown in Equation^ Let [g] denote the linear reformulation of g according to the 
procedure given above. Let us introduce the notation ||_g_|] for the linear lower bound on the polynomial 

lgJ = [g]-E + (g) 

And similarly [j~g]] for the linear upper bound: 

M = [g]-E-(g) 

The construction above guarantees that the linear bounds are correct for all feasible values of the 
original variables x and the corresponding values of the reformulated variables w: 

MM < g(x(w)), Vx(w) g O x (49) 

MM > g(x(w)), Vx(w) g o x (50) 

where x(w) is the point in original variables corresponding to the reformulated point w. The construction 
guarantees that every point x in O x has at least one corresponding point w in O w , and that each point 
w maps to a unique point x(w); however there may be several feasible points w(x) to represent any given 
x. Note that the reformulation [g] is a linear function of the unit variables, remainder variables, and unit- 
product variables which constitute the vector w shown in Equation[5] As each error bound E~(g) and E + (g) 
is a real number, the linear bounds ||_g_|] and \\g~\\ defined above are therefore also linear functions of the 
variables in w. 
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3.3 Pair of Bounding Mixed Binary Linear Programs 

We can use the linear bounds on polynomials described above to generate more and less restrictive versions 
of the original polynomial program PP (inner and outer approximations). As each reformulated program will 
contain several binary unit variables as well as continuous remainder variables and unit-product variables, it 
will be a mixed integer linear program whose integer variables are binary. For the optimistic case in which 
a lower bound on the global minimum is desired, the linear lower bound [[f_|J on the objective function f 
should be used in the optimistic reformulated program LP~. Also, each constraint gj < in the polynomial 
program PP should be replaced in the linear program LP~ with its linear lower bound [[gjJJ ^ 0, which 
is less restrictive. For standardization each constraint gj ^ with the inequality in the opposite direction 
should be replaced with the equivalent constraint — gj ^ 0, and each equality constraint gj = replaced 
with the equivalent pair of constraints gj ^ and — gj ^ 0. Thus the looser problem LP~ is given by: 

minimize ||_f Jj (w) 

subject to |LgiJ|(w)<0, M(w)<0, lLg q J|(w) < 

and w e O w 

Let us use w to designate the point at which the minimum solution to the optimistic program LP~ occurs, 
and []f]J(w~) to denote the value of the linearized version of the objective function at that point. The linear 
reformulation LP~ can be infeasible only if the original polynomial program PP is infeasible. It can happen 
that the relaxed program LP~ is feasible although the original program PP is not. 

Similarly, for the pessimistic case in which an upper bound on the global minimum is desired, the 
objective function and each constraint in PP should be replaced with its linear upper bound to produce the 
tighter reformulated linear program LP J : 

minimize [[f|j (w) 

subject to |rgiH(w)<0, MlwXO, M(wK0 
and w e O w 

Let us use w + to denote the point at which the minimum solution [[f|](w + ) of the pessimistic program LP J 
occurs. The pessimistic reformulation LPj can be infeasible even if the original program PP is feasible; 
however if LPj is infeasible then PP must be infeasible as well. 

We must use care in interpreting the solutions to the reformulated linear programs LP~ and LP J. Let 
us use x* to designate a point at which the global minimum solution f(x*) to PP occurs. Let us say that a 
point w in the reformulated variables is polynomial feasible if the corresponding point x(w) in the original 
variables satisfies the bound constraints oa < Xi < (3| and the polynomial constraints gj ^ in PP. 

The properties of the solutions to the reformulated programs are as follows. The solution to LP~ places 
a lower bound on the true solution to PP (there is no better solution): 

|LfJ|(w-)<f(x*) (51) 

If it happens that w~ is polynomial feasible, then the value of the original polynomial objective function 
f evaluated at the corresponding point x(w~) is an upper bound on the true solution to PP (there is a 
solution at least that good): 

f(x*)<f(x(w-)), ifx(w-) e <5 x andgi(x(w-)) < 0, g 2 (x(w~)) < 0, g q (x(w"))<0 (52) 

However if the solution point w~ to LP^ is not polynomial feasible then the loose reformulation LP~ does 
not provide any upper bound on the true solution to PP; in that case it is necessary to use an alternate 
means to generate an upper bound. One way is to use the pessimistic reformulation LPj to compute an 
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upper bound on the global minimum solution to PP. If the tighter program LP+ is feasible then its solution 
is an upper bound on the global minimum: 

f(x*) < ^(w+) (53) 

If the reformulation LP + is feasible, then by the construction of LP + the point x( w + ) in the original variables 
corresponding to the solution point w + to LP + must be polynomial feasible. Therefore we can use the value 
of the original objective function f at that point as a tighter upper bound on the global minimum: 

f(x*u f(x( W +)K irfn(w + ) (54) 

In the case that LP + is infeasible it does not provide an upper bound on the global solution to PP; and such 
a result does not prove that PP is infeasible. 

Example 7 Let us consider problem PP1 from Example [2] The solution |JfJJ(w~) = —124.799 to the 
optimistic reformulation LPl^ 2 2 i occurs at the point w~ : (uij = 0,111,2 = 1,11^3 = 0,U2j = 0, 112,2 = 
0,11.3^ = 1,U32 = 1,Ti = 0.666667, T2 = 0,T3 = 1). This value [[fj](w~) is a lower bound on the true 
global minimum of PP1. The corresponding point in the original variables x(w _ ) : (xi = 3,X2 = 0,X3 = 8) 
happens to be polynomial feasible, satisfying the bound constraints and two additional constraints in PP1. 
The original objective function evaluated at this point has value f(xi = 3,X2 = 0, X3 —8) — — 119. This 
value f (x(w - )) is an upper bound bound on the true global minimum of PP1. Thus the true global minimum 
solution to PP1 lies in the interval [—124.799, —119]. It turns out that the global minimum to PP1 is indeed 
—119 which we can prove by reformulating the problem using smaller error limits and correspondingly more 
unit and unit-product variables. 

Example 8 Let us consider a different polynomial optimization problem PP2 which is used as Example 2 
in Li and Chang [7 : 

minimize : 0.6224x 3 x 4 + 19.84xfx 3 + 3.1661xfx 4 + 1.7781x 2 x§ 
subject to : xi ^0.0193x3 
x 2 ^ 0.00954x3 

1 .33333x1 7t + x§x 4 7r ^ 750.173 
x 4 < 240 

and: xi £{1,1.0625,1.125,1.1875,1.25,1.3125,1.375} 
x 2 e {0.625, 0.6875, 0.75, 0.81 25, 0.875, 0.9375, 1 } 
47.5 < x 3 < 52.5 
90 sC x 4 < 112 
tt = 3.14159 

Note that we can accommodate the discrete variables xi and X2 by reformulating each of them with binary 
unit variables as in Proposition [l] but without the continuous remainder variable ri or xi. The rest of the 
method works without modification. In this case the reformulation of each variable is given by: 

X! = 1+0.0625x11,1+0.12511.1,2 + 0.2511.1,3 
x 2 = 0.625 + 0.0625u 2 ,i +0.125u 2 ,2 +0.25u 2 , 3 

x 3 = 47.5 + 0.00976562u 3J + 0.0195312u 3 , 2 + 0.0390625u 3 , 3 + 0.078125u 3 , 4 + 0.15625u 3 , 5 

+ 0.3125u 3 , 6 + 0.625u 3 , 7 + 1 .25u 3>8 + 2.5u 3 , 9 + 0.00976562r 3 
x 4 = 90 + 0.6875u 4)1 + 1 .375u 4i2 + 2.75u 4! 3 + 5.5u 4 , 4 + 1 1 114,5 

+ 0.6875r 4 
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using the error tolerances (ki , K2, K3, K4) = (0.0625, 0.00976562, 0.6875, 0) and the corresponding numbers of 
unit variables (o"i , 02, °"3> 04) = (3, 9, 5, 0) The loose reformulation LP2^ 9 5 j of this problem PP2 has the 
solution [|_f|J(w - ) = 6395.51 which is a lower bound on the true minimum f(x*) as shown in Equation [5TI 
This lower bound occurs at the reformulated point w~ : (ui 1 = 0, ui,2 = 0, ui ; 3 = 0,U2,i = 0,U2,2 = 
0,u 2;3 = 0,u 3)1 = 0,u 3 , 2 = 0,u 3 , 3 = 0,u 3)4 = 0,u 3 , 5 = 0,u 3 , 6 = 0,u 3)7 = 0,u 3 , 8 = 0,u 3 , 9 = 0,u 4>1 = 
0, U4 ; 2 = 0, U4 ;3 = 0,114,4 = 0,U4,5 = 0, ri = 0, T2 = 0, r 3 = 0, T4 = 0, Ts = 0). The corresponding 
point x(w _ ) : (xi = 1,X2 = 0.625, x 3 = 47.5, X4 = 90, n = 3.14159) in the original variables happens to 
be polynomial feasible; the original objective function evaluated at this point has value f (xi = 1 , X2 = 
0.625, x 3 = 47.5, x 4 = 90, n — 3.14159) = 6395.51 . Thus the true global minimum solution to PP2 lies in the 
interval [6395.51,6395.51]. This is a better solution than computed in Li and Chang [7J. More importantly, 
the method presented here assures that this is in fact the global minimum as there can exist no better 
solution than the lower bound 6395.51. 

Example 9 We consider the problem PP3 which is Problem 338 in Schittkowski [10] and Example 3 in [7]: 

minimize : — x^ — x\ — x| 
subject to : jX-\ + X2 + x 3 = 1 

A + 1*2 + 4-*3 = 4 

and : — 2 ^ xi ^ 2 

-2.45 < x 2 < 2.45 

-4^X3^4 

(The bounds on the variables were not specified in the original formulation but they are implied by the 
equality constraints and made explicit here.) The program is reformulated using (o"i , o"2, cr 3 ) = (7,7,7) 
and correspondingly (ki,K2,k 3 ) = (0.03125,0.0382813,0.0625). The solution to the loose reformulation 
LP3^ )7)7) is l[fJJ(w - ) = -10.9965 which occurs at x(w") : (xi = -0.375,x 2 = -1.65897,x 3 = 2.84647). 
However this point x(w~ ) is not polynomial feasible; it violates the second constraint in PP3. It happens that 
the tight reformulation LP3jt, 7 7) is infeasible, as is often the case for polynomial programs with nonlinear 
equality (as opposed to inequality) constraints. As an alternate means of finding an upper bound on the 
true global minimum to PP3 (and a polynomial-feasible point at which that upper bound occurs) we can 
generate a focused problem that uses narrower ranges of values for the variables xi , X2, and x 3 concentrated 
near the polynomial-infeasible point x(w _ ) where the solution to LP3(^ 7 7 , occurs. If we choose a range 
focused on x(w _ ) ± K| for each variable then the resulting polynomial program is the following (let us call 
it PP3.2): 

minimize : — x\ — x\ ~ x\ 
subject to : -jXi + X2 + x 3 = 1 

A + hi + \ x l = 4 

and : -0.40625 $$x, < -0.34375 
-1.69725 x 2 s$ -1.62069 
2.78397 < x 3 2.90897 

This program PP3.2 can be reformulated using the same numbers of unit variables (o"i , o"2> 03 ) = (7, 7, 7) used 
for PP3 above, which now produce smaller error limits (ki , K 2 , k 3 ) = (0.000488281 , 0.000598125, 0.000976562) 
due to the narrowed bounds on each variable. The pessimistic reformulated focused program LP3.2^ 7 7 j 
has solution —10.9928 which is achieved at the polynomial-feasible point x(w + ) : (xi = — 0.36621 1,X2 = 
—1.6622, x 3 =2.84531). Combining the results of these reformulations LP3^ 7 7 j and LP3.2p, 7 7 > shows that 
the global solution to the original problem PP3 must lie within the interval [—10.9965,-10.9928]. Again, 
in contrast to the approximation method presented in Li and Chang [7J , the bounding approach presented 
here guarantees that there cannot exist a better minimum than —10.9965. 
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3.4 Alternative Formulation: Allowed Constraint Violation 

Another way to use the reformulation technique described above would be to compute a single mixed binary 
linear program LP from the original program PP, using the linearized version [g] of each constraint g and 
the linearized version [f] of the objective f . The error bounds on each [g] could then be used to calculate 
the possible constraint violation Tt for each constraint gi ^ 0, and from these the maximum possible 
constraint violation t across all constraints could be computed. The interval [z~,z + ] would then contain 
the global optimum z* for the variant of the problem PP in which each constraint is nearly satisfied (within 
the feasibility tolerance x). The user could adjust the cq parameters in the pre-solution phase in order to 
achieve the desired feasibility tolerance. Note that in this alternative version, feasibility or infeasibility of the 
reformulated optimization problem does not guarantee feasibility or infeasibility of the original polynomial 
problem. 

4 DISCUSSION 
4.1 Problem Size 

Let us now consider the number of variables and constraints that must be added to a polynomial optimization 
problem PP in the course of reformulation and linearization as described above. We assume that we have a 
list m = {mi , m.2, . . . , mt) of all the terms used in monomials in the polynomial program PP and we use d 
to denote the greatest degree of any term in m. We use o~ = sup t o~i to denote the largest number of unit 
variables required to reformulate any variable x;.. 

The original program PP shown in Equation [T] contains: 

o n bounded real variables xi , xz, ■ ■ . , x n 

o q polynomial constraints gi ^ 0, gz ^ 0, . . . g q < 0. 

o One polynomial objective function f 

In the reformulated program LP there are: 

o Exactly n remainder variables Ty , T2, . . . , r n which are real variables 
o Exactly o~i + o~2 + • • • + o~ n unit variables which are binary 0-1 variables: 

ui,i u 1>2 ••• u 1)(J , 

U-2,1 U-2,2 ' ' ' U 2,(J2 
U-n.,1 u rt,2 • ' • U-n,(j n 

Let us use cf) to represent the number of binary variables in a reformulated mixed binary linear program 
LP. The text above shows that cf) = YJl=~\ °~i which implies cf> ^ no". 

o Not more than td(o~ + 2) d unit-product variables which are real variables. For each monomial term 

m k e M: 

yk,1,1 yic.1,2 ••• yk,1,n k ,i 
yk,2,1 yk,2,2 ••• yk,2,n k ,2 

yk,s k ,l yk,s k ,2 •■' yk,s k ,n k>Sk 

where each Sk is the number of terms Zkj ,Zk,2j • • • > z k,s k needed to express the k-th monomial term 
mk according to Equation [18] and n.k,t is the number of remainder variables n r in the £-th term Zk,t 
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or one if n r = 0. Note that n r is limited by the maximum degree d of polynomial expressions in PP: 
n r ^d. Note also that the number of terms Sk for each monomial mt is limited by the degree dk of m.k 
and the maximum number cr of binary variables used to represent any original variable: Sk (cr+2) dk , 
using the product in Equation [19] Let us use \\> to represent the number of unit-product variables. As 
shown here \L> ^ td(o" + 2) d . 

o Not more than td 2 (d + 1 )(o~ + 2) d product constraints, as for each product variable y^,t,h added it is 
necessary to add n r (n u + 1 ) constraints to satisfy Proposition [3] and both the numbers n T of remainder 
variables (or one if there are none) and n u of unit variables in any element Zk,e are limited by the 
maximum degree d of polynomials. Let us use p to denote the number of additional constraints. 

o Linearized versions of the original q constraints. 

o The linearized objective function [f]. 

To summarize this using O-notation, the reformulated linear program LP will have a number of variables 
that is 0(tdcr d ) including O(ncr) binary 0-1 variables; LP will also have a number of additional constraints 
(besides the original q in PP) that is 0(td 3 o~ d ). Note that all of the added constraints involve binary 
variables. 

An important consideration is the time required to solve each mixed binary linear program. For a program 
with <J) binary variables, in the worst case a branch-and-bound algorithm would require enumeration of the 
2* distinct combinations of and 1 for each variable, and solving a continuous linear program for each case. 
As ct> is O(ncr) the time required to solve each linear reformulation is 0(2 ncr ) multiplied by the time required 
to solve a standard linear program with n + <£> + il> variables and p ~ 0(td 3 cr d ) constraints. 

4.2 Conclusion 

This paper presented a reformulation and linearization technique to generate an approximate solution to a 
polynomial optimization problem. The approximate solution takes the form of interval bounds on the true 
global optimum. In the reformulation step each variable in the original polynomial problem is replaced by a 
sum of binary and continuous variables, the number of which depends on the error limit specified by the user 
for each original variable. In the linearization step products of continuous variables are replaced by sums 
of those variables, and constraints are added to the problem to limit the differences between those sums 
and products. Bounds on the error introduced by linearization are computed, and with these bounds a pair 
of mixed binary linear programs is created whose solutions bound the solution to the original polynomial 
program. The tightness of the generated interval bounds depends on the error limits specified by the user, 
which also determine the size of each reformulated program and consequently the time required to solve it. 
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